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SUMMARY 

The proof oj the theorem that the elliptical distribution 
of lift over the span is that which will give rise to the 
minimum induced drag has been given in a variety of 
ways, generally speaking too difficult to be readily 
followed by the graduate of the average good technical 
school of the present day. In the form of proof here- 
with, an effort is made to bring the matter more readily 
within the grasp of this class of readers. The steps in 
the proof, briefly outlined, are as follows: 

1. Given a basic distribution of lift across the span 
denoted by (a) with a second supplementary distribution 
denoted by (I). Then it is shown that the induced drag 
of lift (a) in the downwash due to lift (1) and the induced 
drag of lift (1) in the doumwash due to lift (a) are equal, 
ind that the total effect of the small distribution {!) on 
the induced drag will be measured by twice either of these 
small quantities. 

£. Next two small changes are assumed in a basic 
distribution (a). These are represented by (£) and {&), 
and are further assumed to be equal in amount and 
opposite in algebraic sign, thus leaving the original lift 
unchanged in amount, but changed in distribution. 
Under these conditions it is then shown thai in order 
for the distribution (a) to be that for minimum induced 
drag, the change in induced drag due to this small change 
in distribution must be zero. 

S. It is next shown that for any pair of small changes 
such as (1) and (#) the only value of the basic down- 
wash which will meet the condition of step (£) is down- 
wash coTistant across the span. 

4. It is known mathematically that the elliptical dis- 
tribution across the span is that which gives a constant 
value of the downwash and hence as a result of (1), {£), 
(8), tkis must be the distribution which will give the 
minimum value of the induced drag. 

The theorem setting forth the elliptical distribution 
of lift over the span of an airplane as that which will 
give rise to the minimum induced drag is one of funda- 
mental importance in connection with the study of 
many problems in aerodynamic theory. Its proof has 
been given by the use of the calculus of variations and 
otherwise, but hitherto, so far as the present writer 
has noted, in form and procedure adapted rather to 



j the expert than to the graduate of the average good 
technical school of the present day. The present proof 
has been prepared with the view of bringing a better 
understanding of this important theorem more readily 
within the grasp of the technical graduate who, so 
far as mathematics is concerned, has no more than a 
fair grasp of the elements of differential and integral 
calculus. 

We at first state the general aerodynamic elements 
involved in the problem, and with which we assume 
the reader to be familiar. 

notation 

Let x denote any point on the span of the plane, 
reckoned from the center as origin. 
yc£r=lift for a short length of span dx at the point x. 




A P O Q D 

Figc** l 



r=the circulation, or vortex strength at the point x. 
i/?= the induced or downwash velocity at any point 

determined by an abscissa x s . 
6 = the half span. 
V= velocity of plane. 
p= density of medium, 
j D= induced, drag. 

Let the distribution of lift across the span be given 
by any curve as in Figure 1. 

"We make only two assumptions regarding the 
character of this distribution. 

(1) The value of y is zero at each end. 
; (2) The curve ABCD has no abrupt changes in 
I curvature. 

j That is, the value of dyfdx varies continuously along 
I the curve from A to D. 

Then the elementary theory of lift gives y — r V p. 
But V and p are constant for the plane and therefore 
y is everywhere proportional to r and vice versa. 

Consider any point P on the span determined by the 
abscissa x. Then at this point, and generally, 
I dr , 1 dy , 
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is an element of change in r and as such is a measure 
of the strength of the elementary trailing vortex given 
off at such point. 

Now let y=j{x) be the equation to the curve 
ABGD. Then we have 

%dx=f (a) d*and 

Then the theory of induced or downwash velocity 
gives, for the poult Q, 

A 1 dr , 1 f (a) dx 

as the measure of the element of v> at Q due to the 
action of the vortex filament at P. (Reference 1.) 

Then the entire w at Q will be given by the integra- 
tion of the above expression across the span or from 
A to D. This gives: 



_ 1 Cf (; 
W ~4ir VpJ ~ t 



(a) da 



(1) 



Again the theory of induced drag gives, for the ele- 
ment of D at the point Q. 

V) 

dZ? = j/p,dzi 

where y and w both refer to the point Q. Hence: 

D=^>Jywdxi (2) 

Equations (1) and (2) are basic with reference 
to the theorem which we have now to develop and, as 




noted above, it is assumed that the reader is familiar 
with the general features of the vortex theory of lift 
and with the establishment, from fundamental prin- 
ciples, of these two equations. 

We shall now proceed to establish the theorem in 
question by a series of steps which may be listed as 
follows: 

(1) Given a basic distribution of lift across the 
span. Denote this by (a). Given a second small 
or supplementary distribution. Denote this by (1). 
Then we have to show that the induced drag of lift 
(a) in the downwash due to lift (1) and the induced 
drag of lift (1) in the downwash due to lift (a) are 
equal, and that the total effect of the small distribu- 
tion (Jf ) on the induced drag will be measured by twice 
either of those small equal quantities. 

(2) We next suppose two small changes made in a 
basic distribution (a)- These may be represented by 



(1) and (#). These are further assumed to be equal 
in amount and opposite in algebraic sign, thus leaving 
the original lift unchanged in amount, but changed in 
distribution. Then under these conditions we have 
to show that in order for the distribution (a) to be 
that for minimum induced drag, the change in induced 
drag due to this small change in distribution must be 
zero. 

(3) We have next to show that for any pair of small 
changes in distribution such as (1) and (2) the only 




Figure 3 

value of the basic downwash v> a which will meet the 
condition of step (2) is w a = constant, or downwash 
for the basic distribution = constant across the span. 

(4) It is known mathematically that the elliptical 
distribution of lift across the span is that which gives 
a constant value of the downwash w and hence is that 
which will give the minimum value of the induced 
drag. 

We shall now proceed with these steps in order. 

(1) In Figure (2) let (a) represent a basic distribu- 
tion of lift across the span and (1) a second or small 
supplementary distribution. Then in equation (1) 
we may omit the constant factor and write: 



J {xi-x) 



(3) 



Let us integrate this expression in accordance with 
the formula 

fu<Lv=uv—fvdv, 

Where u represents j^hr^) ^ dv tenotesf (x) dx. 
We shall then have — 



where -b and + b denote the limits of the integration 
at A and Z?, Figure 1. 

Taking the first term on the right, special caution 
must be exercised in making the evaluation due to the 
fact that in passing from x ■» — b to x=> + b } the denom- 
inator passes through 0. The range of values of this 
expression will be therefore as indicated in Figure (3). 
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At A the value of /(a) is zero, but a zero approached 
through positive values. We may therefore term 
such a zero a+0 and similarly one approached through 
negative values, a— 0. We shall then have at A, 
J (x){ —%) = +0* Then as the point under considera- 
tion approaches Q from A 9 the value of (x L — x) will 
approach a positive zero and the value of/ fa) (xi—x) 
at Q becomes + ». Then as the point passes beyond 
Q where x>x u the zero shifts from positive to negative 
and the value just on the right of Q becomes — co and 
then gradually decreases numerically as z increases, 
finally approaching 0 on the negative side at D. 

The total change in the value of the expression 
from A to D is therefore from -f 0 to —0 passing 
through + co and — co on the way. Suppose, however, 
that we take two points Q t and Q n respectively on 
the left and right of Q } and at a very small distance 
e from Q. Then if for brevity wo represent the ex- 
pression/Car)/ fa — x) by a single letter tp, we may write — 



Ilm.e-0 



But as we have seen, p A and pd=*0 and hence 



-•I 



tlm.<-0 



iIm.f-0 



But at Qt, f(x)=f(xi-e) and {x t -x)=*e* Likewise 
at Qrf{z) =/(a?i+6) and (pi—x) = — e. Hence we have 
finally, restoring the original form, 

Xi-XJa L * J 

This we may write 1 — 

2/Cnn 

We shall then have for equation (4) 

Now turning to equation (2) we may, for our pres- 
ent purpose, omit the constant V and write 

D~fyw&x x (6) 

Note that in this equation both w and y refer to 
the point 0 and the integration is with reference to a*. 

Applying this equation to the case assumed in 
Figure 2 we note that the total value of D will result 
from four constituent elements as follows: 

(a) That due to lift (a) in downwash (a) 

(b) That due to lift (1) in downwash (1) 

(c) That due to lift (a) in downwash (i) 

(d) That due to lift (1) in downwash (a) 

Then considering (a) as basic, the change will be the 
sum of (b), (c), and (d). But if distribution (1) is 
small in magnitude, item (b), will be small of the 

1 For a somewhat more formal mathematical treatment of tnis evaluation, see 
Appendix 1. 



second order and negligible in comparison with the 
other values. We shall have, therefore, for the re- 
sultant effect of the supplementary distribution (1), 
simply the sum of (c) and (d). 

Now referring to Figure 2, let us denote by sub- 
scripts a and 1 the various quantities referring to 
these two distributions of lift. Then remembering, 
as above noted, that in (6) y refers to the point 0, as 
fixed by the abscissa x u we may write: 

y*=/« (xi) 

and from (5) 

Next we apply (6) taking w a and y u This will give 
the induced dcag for lift distribution (/) in the down- 
wash caused by distribution (a). We thus have 

Hm.e-0 

Or writing the last term as a double integral, 



Viauxa 4 

Then vice versa if we take y a and w t we shall have 
the induced drag for lift distribution (a) in the down- 
wash caused by distribution (1). This will give: 

(9) 

Now in equations (8) and (9) it is' clear that the first 
terms on the right, even though they do involve infini- 
ties, will be equal. For the second terms, we refer to 
Figure 4. The square with center at 0 and axes of x 
and x t represents the combination of values of x and x t 
to be covered by the integration. That is, every pos- 
sible combination of values of x and x t between +6 
and —5 is represented by a point on this square and 
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the double integral means that for each such point or 
combination, the value of the expression under the sign 
of integration is to be found, and these elements 
summed. Now at a point P we shall have, in equation 
$)fafa) =MAP) and/, (z) =MOA), while at the point 
Q we shall have, in equation (9), jtix^^j^OB) and 
fa(?) = fa(BQ). But for any point P there will always 
be another point Q for which AP=BQ and OA=OB. 
And what is true in the first quadrant will be equally 
true in the others. Hence for every point represent- 
ing a value in the numerator of equation (8) there will 
be another point representing a value in equation (9) 
such that the two values are equal, each to each. The 
denominators likewise will be equal since {x t — x) will 
have the same numerical value for each point, and in 
the square the value will be the same regardless of the 
sign. Hence the summations will be equal and the 
values of (8) and (9) are equal throughout and as a 
result, we have, Dai=D la or — 

The induced drag due to y t under downwash w a will 
equal that due to y a under downwash vi u And further- 
more since, as we have seen, the total effect is measured 
by the sum of these two items, it will obviously be 
measured by twice either one taken singly. 

It may be noted at this point that the proof that 
Dia^Bai is independent of the relative size of distribu- 




2 

Figure 5 



tions (a) and (1) and that this relation is true generally. 
For our present purposes, however, we are concerned 
rather with the particular case where (1) is small and 
where the induced drag due to yt in the field w f is 
negligible and .\ where the total effect is measured as 
above noted. 

(2) The second stage in the proof calls for a con- 
sideration of the results of making two small changes 
in a basic distribution of lift, one positive and one 
negative, equal in numerical value and of any char- 
acter of distribution: See Figure 5. Here (1) and (#) 
denote any two small distributions of lift, No. (1) 
positive and No. (2) negative, equal in numerical value 
and distributed each in any manner along the span. 

Due to the fact that (1) and (#) are both small it 
will result that the elements of induced drag due to 
y s v) t , yiV) t , y»v)t and y t w g will all be small of the second 
order and therefore negligible in comparison with the 
other values. It will result then that the total change 
in induced drag will be measured by the components: 

D at = drag due to y a in field of w t 
D Ia « drag due to y t in field of w t 
Dfl^drag due to y a in field of w* 
Z?* a = drag due to y t in field of w a 



But as we have seen, 

Dai = Di a and Dat=D ia and hence 
+ D Ia +1)0,+ Z>„ = 2{D Ia + D ia ) 

That is, the total change in the induced drag as a result 
of these changes in the distribution of lift (the total 
lift remaining the same) will be twice the sum of the 
induced drags due to y t and y t} both in the down- 
wash ?0 a . 

Now this total change may be (— ), (+), or 0. If 
it is (— ), it will show that the original distribution 
y a was one which admitted of reduction in the induced 
drag for the same total lift, and hence that it was not 
the distribution for minimum induced drag. 

If it is (+), it would be possible, by interchanging 
the signs of the distributions (1) and (£), to make it 
(— ); that is, by making distribution (1) negative and 
(#) positive, we should interchange the numerical rela- 
tions and the result would change sign, This again 
would show that the original distribution y m is not that 
for the minimum induced drag. 

If then the value 2(D la + Dg a ) is anything other than 
zero, it will indicate that the basic distribution y a is 
one which will give an induced drag permitting of 
reduction (by redistribution) by an amount at least 
measured by this value. Hence the only distribution 
of y a which will give an induced drag not permitting 
of reduction is that which will give zero for the value 
of the change 2(JD ta +D $a ). But the condition for this 
is —D ta . From equation (6) this will mean: 

fyiWa da, « - Sy%v) a Ax t (10) 

while at the same time we have the condition: 

fVt Ax^-fytdxt (11) 

(3) To carry out step (3) of the proof, we note that 
the simultaneous fulfillment of equations (10) and (11) 
may be realized in two ways. 

(1) For any and all pairs of distributions y t and y t 
which meet equation (11) the only value of the basic 
downwash w* which will meet equation (1 0) is 
w a = constant. 

(2) For any particular pair of distributions y t and y, 
fulfilling equation (11), distributions of w* may be 
found, not necessarily uniform in value, meeting equa- 
tion (10) for this special pair y i} y, but failing for other 
pairs which still meet equation (11). 

In order not to interrupt the main argument at this 
point, "further reference to this point is transferred to 
an appendix. 

It appears, however, that the fulfillment of equation 
(10) for any and all possible pairs of distributions y t 
and y» which meet equation (11) (as is implied by the 
character of the preceding argument) can only be 
realized if w a is a constant multiplier under the sign 
of integration. 
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That is, in order that the original distribution y € 
shall be that for TrnnrrrniTn induced drag, it must be 
such a distribution as will give a constant value of the 
downwash to across the span. Under these conditions 

then, a small supplementary induced drag ^fyiwAxt, 

will equal numerically a second small supplementary 

induced drag ^fy*w&Xt, provided J % yAx 1 = J % y t dx 1 ; 

and as a result, where y t and y t are opposite in sign, 
wo shall have D ttL = —T>u or Dm+De^Q, or no change 
in the over-all induced drag for this small change in 
distribution. This condition, as we have seen, implies 
that the basic induced drag is at its minimum value 
for the total lift as given by the distribution y a . 

(4) Now it is known that the elliptical distribution 
across the span is that which will give a uniform value 
of the downwash velpcity That is, it is known that 

for the elliptical distribution the value of f ^ fe)^ 

is a constant and equal to x. It is no part of the pres- 
ent paper to prove this particular result. It is a mat- 
ter of integral calculus, pure and simple, and while the 
integration requires special care due to the fact that 
in passing from — b to +fi the denominator becomes 0 
and infinite values are involved, it is here assumed that 
this particular result is either known or accepted. 
(Reference 2.) 

The proof is thus complete that the ellipitcal distri- 
bution of lift across the span is that which will give 
the minimum value of the induced drag for the given 
total lift so distributed. 

An interesting side light on this problem may be 
obtained by considering the results of a continued 




total distribution from something like that of (a) over 
into (6) (See Fig. 6) and passing through the elliptical 
form at a mid-point as indicated by the dotted line. 
Such a series of changes could, of course, be realized 
by the continued application of a series of equal plus 
and minus small changes. In Figure 7 let a a; denote 



X 



series of small changes in distribution (the total re- 
maining the same) such as to gradually change the 



0 E A 

FlOUBK 7 

an axis of change and let A be the point for distribu- 
tion a, B the point for distribution b and 0 that for the 
elliptical distribution. Then AP may represent the 
value of the induced drag for distribution a, BQ for 
distribution 6 and OB for the semi-ellipse. Then the 
entire series of changes by which a is transformed into 
5 will give a series of values of the induced drag as indi- 
cated by the curve PBQ and it is seen how any small 
change from a distribution represented at E will 
result in a slight decrease or increase in the value of 
the ordinate EF depending on whether the change is 
toward or away from the elliptical distribution. This 
will indicate, as we have seen, that EF is not a mini- 
mum ordinate and that the distribution of lift indi- 
cated at E\s not that for minimum induced drag. On 
the other hand, with the distribution at 0 (elliptical) 
and the value of the induced drag 0B } any small 
change either way will make no sensible change in the 
ordinate OB; or otherwise the sum total of the change 
in OB will be zero and this will indicate that OB is a 
minimum ordinate and that the distribution at 0 
(elliptical) is that for the minimum value of the in- 
duced drag. 
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For the evaluation of the expression f(p)/{xi—x) 
between the limits of —6 and +6, we note first that 
the function is assumed to be continuous and without 
singularities. At any point near Q, distant e, we 
may therefore develop f{x)=*f(xi±e) by Taylor's 
theorem. 

Thus for the point Q r on the right of Q we shall 
have: 

m -/W +f(x 1 )e+£^e*+etc. 

And for the point Q x on the left of Q we shall have: 

/(*) -/' e? - etc. 

Hence denoting by <p the function to be evaluated, 
we shall have: 



1 . /fa) f{x,)e f'fa) 
J Or Xi — x Xi—x 2 fa — a) 

Jgi Xi—x xi-x 2(xi-x) 

But for Q TS {xi—x) = —e and for Qfc fri— ar) = +e. 
Hence we shall have: 

Toward the limit the terms having e in the'numerator 
with powers of 1 and higher will vanish and we shall 
have: 

J - J 1 - 2 /fa)l 

&^ QJ <P Jorjllm- «-0 — — ~ J 1Im " e ~° 
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Conditions under which we may have — 

fy L dx = fy 2 dx and (1) 

fyiwdz^fyzwdx (2) 

In Figure 8 let y x and y 2 fulfill equation (1). Then 
such condition may be expressed by 




Area ABC- Area GDE 

Now it is obvious that as we multiply the ordinates of 
y t and y* by w as a factor, we shall likewise multiply 
the ordinates of these difference areas by the same 
factor w. Then in order to insure the fulfillment of 
equation (2), we have only to insure the multiplication 
of these areas ABC and CDE, as a whole, by the same 
factor. Thus if we multiply the areas ABC and GDE 
both by some number m, the enlarged areas will still 
be equal numerically, and equation (2) will be fulfilled. 
Now we can multiply these areas by m by multiplying 
each'and every ordinate by m, or otherwise by various j 
41630—81—25 



other arrangements of the factory such however that, 
as a whole, the areas ABC and CDE will both be multi- 
plied by the same over-all factor m. These special 
distributions of w all depend on the particular pair of 
distributions y u y 2} and would not hold in general for 
other distributions which still fulfill equation (1). 

On the other hand if w is a constant factor, then no 
matter what the pair of distributions y t y 2j equation 
2) will be fulfilled. In such case we can, of course, 
take the w out from the sign of integration and write 
for equation (2) 

wfyi&x = wfy-idx 

which is the same as (1) 

But it is clear from the manner in which these small 
supplementary distributions y t and y t are used in 
developing the line of proof, that they are entirely 
unrestricted as to form and character, so long as 
equation (1) is fulfilled. Or otherwise, to any one 
basic distribution giving a downwash denoted in 
general by id, we may, without in any way affecting the 
argument, apply any and all types and forms of small 
supplementary distribution yi and y 2 , so long as 
equation (1) is fulfilled, and the only single distribution 
of w which will satisfy equation (2) with any and all 
of these possible pairs y t and y 2 , is w= constant. And 
thus this link in the general proof is established. 
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For those who may be interested, the following brief 
outline of the formula, lift = TpV is given. The de- 
tails may be found in any good text on fluid mechanics. 

Given a circular cylinder of infinite or indefinite 
length placed in an infinite stream of fluid having two 
component motions: 

(1) A horizontal motion with velocity 7 in a 
direction _L to the axis of the cylinder. Assume this 
to be from right to left. 

(2) A motion of rotation (assumed left handed) 
about the axis of the cylinder with velocity at any dis- 
tance r measured by F/2xr. 

Then for any circle with radius r, the product of the 
length of the circumference by the velocity will be the 
constant r. This is called the circulation, or vortex 
strength per unit length of the cylinder. 

For an element of length di, the corresponding 
vortex strength will be TcLc. 

Let p be the density of the fluid. Then for the total 
force reaction on the cylinder per unit length, the 
principles of the mechanics of fluids (see any text- 
book) give: 

F=T P V 

The theorem shows also that this force is directed JL 
to the direction of flow of velocity V and such that the 
body is urged toward the side on which the velocity 
due to the circulation, or vortex motion, tends to 
augment the direct velocity V. With the above as- 
sumptions as to directions, this will be vertically up- 
ward. Such a force is called a "lift." 

We next assume a geometrical or imaginary cylin- 
drical surface of radius R surrounding the given 
cylinder of radius r. The total reaction between the 
fluid outside this surface and that inside will consist of 
two parts: 

(1) A net resultant pressure acting over this surface 
and directed the same as the net force on the cylinder 
itself. 

(2) A change of momentum directed downward ex- 
perienced by the fluid during its passage through the 
volume lying between the two cylindrical surfaces. 

But the time rate of the change of momentum is the 
measure of a force and the reaction of this change of 
momentum downward will be a force upward to be 
considered as acting over the surface of the cylinder of 
radius R. 

The same principles of mechanics will then show that 
the sum of these two force reactions is again TpV. The 
two parts have values as follows: 



where F p denotes the part due to pressure direct and 
F n that due to change of momentum. 

If R=r 7 F m =0 and the entire reaction is direct 
pressure as we have seen. If R is very great, the two 
parts approach equality at the half value TpVf2 and 
at co these values are reached. The sum, however, as 
shown by the values in (1) and (2) is always TpV. 

Now, suppose R very large in comparision with r. 
Then imagine the cylinder with radius r to be quickly 
withdrawn and an indefinitely long body with airfoil 
section substituted in its place, the latter being of 
the same general order of size as the cylinder of radius 
r. Or otherwise, consider a separate case, with the 
airfoil section instead of the cylinder. With R very 
large compared with the dimensions of the body, it 
can make no sensible difference in the character of the 
flow at and through the surface of the outer cylinder. 
Or otherwise, we may say that we can certainly go to a 
radius R so great that the flow at and through the 
surface of a cylinder of this radius will differ insensibly 
in the one case as compared with the other. 

But over this surface of radius R, the total reaction 
(pressure + momentum change) will be TpV. 

Consider again the condition of the fluid between 
the two cylindrical surfaces. It is constantly under- 
going a change of momentum as indicated in equation 
(2). Furthermore, while the force reaction of the 
fluid which is thus undergoing a change of momentum 
is upward,- or in the same direction as the force F f} the 
direction of the change of momentum itself is down- 
ward, or opposite to the direction of F 9 . 

Now this fluid can be acted on from the outside, 
only by way of pressure transmitted to it over two 
surfaces — , the surface of the body of airfoil section 
and the surface of the cylinder of radius R. The 
final downward momentum must then be the result 
of the joint action of these two pressures. Call the 
total force reaction, body to fluid, F. Since the 
pressure F 9 is up, the force F must be down and we 
shall have 

F-F 9 =F n or F=F v +F M =TpV 

This gives the total force, body to fluid, down and 
hence the equal force, fluid to body, is upward, and 
measured by TpV, the same value as for the body 
with circular cross section. Hence regardless of the 
form of cross section of the body, the total force 
reaction under the conditions assumed will be a lift 
measured by TpV, 
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The effect of downwash is to produce a vertically 
downward current of velocity w about the airplane, 
compounding with the relative horizontal velocity F. 
The result is a direction of resultant Sow relative to 
the plane, tipped down at an angle, whose tangent is 
wfV. Denote this angle by y. But the lift as defined 
and as given by formulas and by measurement is al- 
ways at right angles to the direction of relative wind. 
Hence with downwash, the direction of the lift L will 
be tipped back from the vertical at the same angle 7. 
But lift in this direction will have a backward compo- 
nent L sin 7. The angle 7 is always small, however, 



and we may, therefore, without sensible error, take the 
tangent for the sine and write L sin y=Lw[V* This 
force opposing the motion of the plane, is called the 
induced drag- For a small element dz, of the span, it 
will have the value 

ID 

<W=Uy&Xi as in equation (2) 



Stanford University, 

October 12* 1929. 
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